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ABSTRACT 

Context. The present work deals with the detection of phase changes in an exoplanetary system. HD 46375 is a solar analog known 
to host a non-transiting Saturn-mass exoplanet with a 3.0236 day period. It was observed by the CoRoT satellite for 34 days during 
the fall of 2008. 

Aims. We attempt to identify at optical wavelengths, the changing phases of the planet as it orbits its star. We then try to improve the 
star model by means of a seismic analysis of the same light curve and the use of ground-based spectropolarimetric observations. 
Methods. The data analysis relies on the Fourier spectrum and the folding of the time series. 

Results. We find evidence of a sinusoidal signal compatible in terms of both amplitude and phase with light reflected by the planet. 
Its relative amplitude is AF p /F+ = [13.0,26.8] ppm, implying an albedo A = [0.16,0.33] or a dayside visible brightness temperature 
7b - [1880,2030] K by assuming a radius R = l,lf?j up and an inclination i = 45°. Its orbital phase differs from that of the radial- 
velocity signal by at most 2o" RV . However, the tiny planetary signal is strongly blended by another signal, which we attribute to a 
telluric signal with a 1 day period. We show that this signal is suppressed, but not eliminated, when using the time series for HD 46179 
from the same CoRoT run as a reference. 

Conclusions. This detection of reflected light from a non-transiting planet should be confirmable with a longer CoRoT observation of 
the same field. In any case, it demonstrates that non-transiting planets can be characterized using ultra-precise photometric lightcurves 
with present-day observations by CoRoT and Kepler. The combined detection of solar-type oscillations on the same targets (Gaulme 
et al. 2010a) highlights the overlap between exoplanetary science and asteroseismology and shows the high potential of a mission 
such as Plato. 

Key words. Stars: planetary systems, Techniques: photometric, Methods: data analysis, Stars: individual: HD 46375 



!> ■ 1 . Introduction lightcurve and detect or place uppe r Umits on the brightness of 

the planet (Guillot & Vannier 20QJ- 
The star HD 46375 is a relatively bright (m v = 7.89, M bol = 1, 
H ! 5.04) K0 star, known to host a Saturn-mass planet (M p sin i = L The . detection of reflected light from an exoplanetary atmo- 
f*? ■ 0.226(19) M Jup ) orbiting in 3.023573(65) days, at 0.0398(23) s P here 18 ver y lnterestln g for several reasons: the amplitude of 
AU, on a quasi-circular orbit e = 0.063(26) dMarcv et al.ll2000l the modulation helps us to constrain the atmospheric proper- 
Bud eTitffl 120 06). The planet does not transit its host star tles ' and ^ particular the planetary albedo. This has direct con- 
(iHenrvl 120001) . HP 46375 is the brightest star with a known sequences for the calculation of the thermal evolution of these 
close-in planet in the CoRoT accessible field of view. As such, P lanets and therefore their structure. Departures from a sinu- 
it was targeted by the CoRoT additional program and observed soidal shape may be caused by non-isotropic scattering and be- 
in a CCD normally de dicated to the asteroseismology pro- tray the presence of small particles in the atmospheres (e.g., 
gram (Baglin et al. 2002), to obtain an ultra-precise photometric I Seager et al. || 2000|) . Finally, the asymmetry of the modulation 

* or a phase shift may be due to horizontal variations in the cloud 

coverage and/or atmos pheric temperatures and trace the atmo- 

Send offprint requests to: P. Gaulme heric circulation ( Showman & Guillot 2002). 

The CoRoT space mission, launched on 2006 December 27, was 

developed and is operated by the CNES, with participation of the In the optical domain, this detection is extremely challeng- 

Science Programs of ESA, ESAs RSSD, Austria, Belgium, Brazil, ing because of the faintness of the signal, of the order of 10 4 

Germany and Spain. ( 1 00 ppm) in the most favorable cases and was achieved in only a 
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few cases. Sing & Lopez-Morales (2009) detected the secondary 
transit of OGLE-TR-56b in the z'-band, a case for which the 
planet is so hot that thermal emission from the plan et contributes 
to the 363 ±91 ppm signal. Using CoRoT data. ISnellen et al] 
(2009) identified the changing phases of the planet CoRoT-lb 
(with an amplitude ~ 130 p pm) and its s e condary transit was 
independently measured bv lAlonso et al.l (l2009bl) . Secondary 
transits of CoRoT-2 b (a 60 ± 20 ppm sign al) were dete c ted by 
lAlonso et all d2009al) . With Kepler data, iBorucki et ail ([2009) 
measured the secondary eclipse of HAT-P-7b (130 ± 11 ppm). 
With HD 46375, we have the opportunity to the detect plane- 
tary phase changes of a non-transiting exoplanet. Such a kind 
detection was achiev ed in the infrared with the Spitzer Space 
Telescope for v And (Harrin gton et al.l l2006) (amplitude ~ 2900 
ppm) but it would be the first time for optical wavelengths. 

In a separate paper (Gaulme et al. 2010a, A&A, submit- 
ted), we used spectroscopic data and seismic analysis of the 
CoRoT lightcurves to constrain the fundamental parameters of 
HD 46375. The star may be considered as a solar analog, of mass 
0.97+0.05 M Q , age 2.6+0.8 Gyr, metallicity [Fe/H] = 0.39+0.06 
dex, and temperature T e g = 5300 ± 60 K. The stellar rotation is 
expected to be in the range [37, 51] days. The spectrometric and 
asteroseismic results provide an estimate of the inclination of 
the stellar rotation axis i m - 50+18° and the planetary mass 
MpSin/ = 0.234 ± 0.008 M Jup . By assum ing that planet orbit s 
in the stellar equatorial plane (see howeve nHebrard et alJ uOlO), 
yields a true planetary mass M p = OJO^jH Mj up . 

In this paper, we use CoRoT observations to search for 
light from the planetary companion. We present successively the 
expected performance (Sect. 2), then, the data analysis in the 
Fourier and the time domain to ensure an ambiguous detection 
of the planetary light (Sect. 3). In Sect. 4, we identify a spuri- 
ous systematic signal in the CoRoT light curve and propose an 
approximate correction. We conclude that we achieve a possible 
detection of the planetary contribution to the light curve, which 
would need further and longer observations to strengthen the re- 
sult (Sect. 6). 



2. Expected signal 

We neglected the eccentricity of the orbit and assumed that the 
starlight was reflected isotropically over the "dayside" of the 
planet (Lambert sphere). The peak-to-peak amplitude of the pho- 
tometric fluctuations caused by the planet is 



AF p (0 _ 2A (R 
F*(t) 
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where F v and F*(f) are the planetary and stellar fluxes, A and 
R the planetary Bond albedo a nd radius, a the o rbital distance, 
and i the orbital inclination (e.g. Charbonneau et alj!999h . There 
are therefore three degenerate parameters: the albedo, the orbital 
inclination, and the planetary radius. 

The planetary radius may be estimated from evolution calcu- 
lations that are param eterized to rep roduce the measured radii of 
transiting exoplanets (lGuillotll2008l) . The result depends mostly 
on the unknown core mass: for an assumed planetary mass of 
0.3 Mj up and a dissipation of 1 % of the incoming stellar light 
at the planet's center, we derive planetary radii ranging from 
0.8 Rj U p for a 60 M ffi dense core to 1 .6 Rj up for a planet with no 
core. 

By assuming that the orbital plane is tilted by i = 45°, 
we found the radius R p = 1.1 R j up and given the large as- 
sumed uncertainty in the planetary albedo, A = 0.05 to 0.5, 
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Fig. 1. Power density spectrum calculated with DFT, with an 
oversampling factor of 8. The dot-dashed curve indicates the 
mean stellar activity. The grey region indicates the range in 
which we expect to detect the planetary signal. 



the relative amplitude AF p /F + of the expected signal is in the 
range [0.8, 61] ppm. Such relative fluctuations of intensity have 
to be compared with the photon noise. In a single channel, 
the amount of photoelectrons collected by the CCD is about 
N* = 1.3 X 10 9 x l0-° An, y s -1 . Hence, for a m v = 7.9 star, the 
relative photon noise is 1/ V^V* = 10~ 3 per second. If individ- 
ual measurements are binned in 6 samples per orbital period, the 
noise level reaches 1 .5 ppm for each orbital sample, for a time 
coverage of 1 1 orbits (i.e., 33 days). Hence, the expected signal- 
to-noise ratio of the planetary signal is [0.5,40] (peak-to-peak). 

Therefore, in terms of photon noise, the precision of CoRoT 
should be sufficient, except for a combination of the lower 
bounds of both inclination and albedo. Note that we have not 
investigated the case w here the reflection is anisotropic (e.g. 
iGuillot & Vannierl l2003). A hypothetical non-strictly sinusoidal 
planetary signal would include harmonics of the orbital fre- 
quency. Since their values are a priori unknown, and given that 
the stellar signal itself displays components of relatively high 
and variable amplitudes at frequencies higher than the orbital 
one, we estimate that the planetary effect will only be detectable 
in terms of its main orbital frequency (Eq. 1). 



3. Ambiguous result from data analyses 

3.1. No clear planetary signature in the Fourier spectrum 

The CoRoT light curve (Fig. ??) does not exhibit signs of high 
activity. Only relative variations smaller than 400 ppm are ob- 
served on typical timescales of 10 days. On the one hand, the 
lack of activity aids the identification of the planetary signal em- 
bedded in the lightcurve. On the other hand, the low activity cou- 
pled with the long stellar rotation period ensures that it is impos- 
sible for the spot an alysis to be performed in a similar way to 
Moss er et alj (12009). Therefore the stellar rotation estimate can- 
not be refined, complicating the disentangling of the planetary 
from the stellar signal. 

We first search for signatures of the planetary chan ging 
phases in the Fourier spectrum. From iButler et"aT] (j2006), the 
planetary orbital period is 3.023573 ± 0.000065 day, which cor- 
responds to the frequency 0.330735 ± 0.000007 cycles per day 
(c/d). The frequency resolution of a 34-day run is 0.029 c/d. 
Hence, the planetary signal should appear as a peak in the fre- 
quency range [0.32, 0.35] c/d. In Fig. Q~l we present the power 
density spectrum calculated using the discrete Fourier trans- 
form (DFT). The background, whose amplitude is the result of 
the stellar activity, was reproduced by fitting a sum of 2 semi- 
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Fig. 2. Folded time series over 3.023573 days. The 2 red sine 
curves represent the minimum and maximum limits of the ex- 
pected signal according to the ephemeris of Butler et al. 2006. 
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Fig. 3. Folding of the HD 46375 and HD 46179 time series over 
3 days. The grey zones correspond to region where the correla- 
tion between both light curves is maximum. 



Lorentzian functions to the whole power spectrum (Gaulme et 
al. 2010a). 

No excess power appears in the expected frequency range. 
Moreover, we note the peak at 0.38 c/d, whose origin is not 
explained by typical systematic signals of the CoRoT, such as 
the orbital period (13.97 c/d), the Earth rotation period (visible 
around 1 c/d), or a combination of them. Therefore, the Fourier 
analysis is inconclusive and a direct analysis of the time se- 
ries has to be performed, by taking into account the information 
about the planetary orbital phase. 

3.2. A sine trend in the folded time series 

The alternative way to search for the planetary signal, for which 
both phase and period are known, consists of folding the time 
series over the planetary period. To increase the signal-to-noise 
ratio, we first rebinned the data over 2 CoRoT orbital periods, i.e. 
6184 s. Then, the low frequency trend, consisting of the moving 
average weighted by a triangular window of 7-day FWHM, is 
subtracted from the time series. We note that we checked this 
step does not affect the planetary period. The rebinned and flat- 
tened time series is folded over the planetary period. If the noise 
is dominated by photon noise, the standard deviation of the time 
series is 176 ppm (Fig. ??) for 92 140 data measurements gives 
for the folded time series (43 bins) has a 3.9 ppm noise level. 

The analysis of the folded time series over the planetary orbit 
detects two main trends (Fig.0. Firstly, a sine trend with a peak- 
to-peak amplitude of about 20 ppm appears to be in phase with 
the expected signal. Secondly, a sharp excess power is visible in 
the range t = [2.2, 2.5] days, which has an amplitude of about 35 
ppm. We must explain the origin of the tall sharp peak, to before 
we can be certain that the sine trend is not an artefact. 
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Fig. 4. Top: folding of HD 46375 time series over the plane- 
tary period. Middle: folding of HD 46179 over the same period. 
Bottom: folding of the time series made of the difference be- 
tween HD 46375 and HD 46179 data. The red sine curves in- 
dicate the data fitting by excluding the regions indicated in Fig. 
® 

4. Evidence of a signal compatible with the planet 

4.1. A 24-h periodic systematic signal 

The peculiarity and difficulty of HD 46375b is that its orbital pe- 
riod is close to an integer number of days, which makes its detec- 
tion very sensitive to telluric systematic signals. Contaminating 
signals with 1-day periods wer e shown to be pre sent in the 
CoRoT exoplanetary lightcurves (Mazeh et al. 2 0091) . 

In our case, among the 9 stars followed up by CoRoT in 
the same observation run in the seismic field, only HD 46179 
('"Hip = 6.68) exhibits less activity than HD 46375, with varia- 
tions below 165 ppm after removing a linear trend. We processed 
the light curve in the same way as HD 46375, and then folded 
both time series exactly over 3 days (Fig[3j. The 2 curves appear 
to be strongly correlated with a clear 24-h period. Therefore, a 
common systematic signal is present in both lightcurves. This 
feature is probably related to the CoRoT orbit, and may be 
caused by e.g. flux variation s generated by the e ntrance/exit in 
the night side of the Earth dAuvergne et al.ll2~009l) . We note that 
only HD 46375 exhibits a sine trend when folded about the plan- 
etary period, which indicated that the sine trend is characteristic 
of HD 46375 (Fig.©. 

4.2. Cleaning of the time series 

We first tried to clean the lightcurves using the SYS -REM algo- 
rithm (Guterman & Mazeh, private comm.). However, this did 
not lead to a reduction in the standard deviation at low frequen- 
cies. We therefore adopted a simple cleaning method, by con- 
sidering the HD 46179 data as reference. Firstly, we flattened 
the HD 46179 lightcurve by removing a linear trend, and then 
subtracted it from the HD 46375 time series. The rebinned and 
folded time series exhibits a clearer sine trend, while the relative 
amplitude with respect to the sine amplitude is reduced by about 
2 (Fig.©. 
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Fig. 5. Power density spectra of the time series corresponding 
to the difference between HD 46375 and HD 46179 lightcurves, 
compared to HD 46375. 

Moreover, the Fourier analysis confirms that while the spec- 
trum amplitude decreases at both 1 c/d and 0.5 c/d, it increases 
significantly in the frequency range which including both the ter- 
restrial harmonics (1/3 c/d) and HD 46375b orbital period (Fig. 
[5]). However, since the light due to Earth-shine is unlikely to fol- 
low a simple path to reach the detector, we cannot reject the hy- 
pothesis that the increase in power at 0.33 c/d could be an arte- 
fact, introduced by the comparison of two light curves obtained 
on different parts of the detector. 

4.3. A signal compatible with HD 46375b 

The variability in the 1-day period spurious signal in both time 
and space prevents us from properly modeling it and cleaning the 
data. From Fig. [3] the intensity of this spurious signal reaches a 
maximum in 0.3-day long intervals close to phases [0.3, 1 .3, 2.3] 
days. We therefore fitted the folded time series by excluding 
these 3 ranges, to minimise the influence of the systematics on 
the sine fitting. The fitting function is a simple sine, whose phase 
</>, amplitude AF p /Fi,, and offset are free. The harmonics of the 
planetary orbit are left aside, because of the low signal-to-noise 
level. 

We first fitted the raw, then the "roughly cleaned" light 
curves. The least squares fitting (reduced^- 2 = 1.75) measured 
a relative peak-to-valley amplitude AF V IF* = 13.0 + 5.0 ppm 
and a phase shift of 19 ± 9° with respect to the e phemeris. A 
phase shift of 19° corresponds to 1.4-cr with the iButler et al.1 
(2006) error bar on the elongation date, which is reasonable 
since their value lies on 50 measurements, all performed before 
2006. The fitting of the "cleaned" data allows us to improve 
the fitting (y 1 = 0.87) and to re-estimate the sine amplitude 
AF P /F* = 26.8 + 3.0 ppm, with a phase shift (f> = 29 + 6°, i.e. 
corresponding to a 2-cr error relative to the predicted ephemeris. 
Moreover, the amplitude of the peak corresponding to the plan- 
etary period is significantly larger in the density spectrum of the 
cleaned data (Fig. [5}- 

After fitting the folded time series by excluding three 
regions, we identified a signal compatible with the planet. 
Assuming that i = 45°, R = l.lRj up , and AF p /F* = 
[13.0,26.8] ppm, our inferred albedo is between 0.16 and 
0.33, which is higher than expected on the basis of theoretical 
models, but within th e range of actually derived albedos (see 
ICowan & Agolll2()10h . Slightly larger radii and inclinations, of 
course, yiel d lower albedos cons istent with theoretical expec- 
tations (see ISudarskv et all 120031) . Alternatively, the observed 
flux may be produced by direct thermal emission. The bright- 
ness temperatures can be calculated from the ratio of the black 
body emissions of the planet with respect to the star, weighted 



by the CoRoT transmission function (lAuvergne et alJl2~009h . We 
found that the minimum temperature is at least [1880,2030] K 
for AFp/F* = [13.0, 26.8] ppm, by assuming that the night side 
has a negligible emission and the planet is seen almost edge- 
on. This is unlikely, because the ratio of dayside optical bright- 
ness temperature to equilibrium temperature would be at least 
between 1 .6 and 1 .8 for HD 46375b compared to measured val- 
ues of between 1.2 and 1.45 for OGLE-TR-56b, CoRoT-lb, and 
CoRoT-2b. 

5. Discussion and conclusion 

Our analysis of observation data presented here has been the 
first attempt to directly detect the stellar light reflected by a 
non-transiting planet in the visible domain. In parallel to the 
present work, we have performed a seismic analysis on the same 
data set. Coupled with ground-based spectropolarimetric obser- 
vations, obtained with the NARVAL spectropolarimeter at the 
Pic du Midi Observatory, we have been able to strongly refine 
the stellar modelling (Gaulme et al. 2010a). 

We used two approaches to detect the planetary contribution 
to the CoRoT light curve. On the one hand, the orbital signal was 
only barely detected by our Fourier analysis. On the other hand, 
the time series analysis, consisting of rebinning and folding the 
light curve over the planetary orbital period, has allowed us to 
detect a signal compatible with the planetary orbital signature. 
The folded light curve appears to be composed of two mains 
signals: a sine trend, in phase with ephemeris, and a high excess 
of intensity, which has no a priori explanation. 

We compared our data with the time series of HD 46179, 
which is the only star among those of the same seismic CoRoT 
run to exhibit a low activity, which is actually lower than that 
of HD 46375. The 3-day folding of both time series enabled us 
to detect a common 1-day period systematic signal, with =s 35 
ppm maximum amplitude (Fig. [3). We reduced by a factor 2 the 
relative amplitude of the spurious peak with respect to the sine 
trend, by subtracting the filtered time series of HD 46179 from 
the HD 46375 time series. 

We illustrated that the data is affected by the presence of a 
systematic signal presumedly of telluric origin, which is strongly 
variable in both time and space. It appears necessary to study 
the photometric fluctuations at low frequency in the whole 3- 
year CoRoT set data especially in the seismic field, to refine the 
modeling of spurious signals. This objective goes far beyond the 
aims of the present work. However, the only way of determin- 
ing whether that the sine trend is of planetary origin consists of 
observing the HD 46375 system for a duration long enough to 
disentangle the planetary orbital period from 3 days, i.e. a dura- 
tion longer than ~ 100 days, as achievable with a CoRoT long 
run. 
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